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SUBARU WEAK-LENSING STUDY OF A2163: BIMODAL MASS STRUCTURE 
N. Okabe\ H. Bourdin^, P. Mazzotta^'*, AND S. Maurogordato'' . 



o 

(N 

CD 



O 
U 

6 



(N 
> 

O 
O 

o 

l> 

o 



% 



ABSTRACT 

We present a weak-lensing analysis of the merging cluster A2163 using Subaru/Suprinie-Cam and 
CFHT/Mega-Cam data and discuss the dynamics of this cluster merger, based on complementary 
weak-lensing, X-ray, and optical spectroscopic data sets. From two-dimensional multi-component 
weak-lensing analysis, we reveal that the cluster mass distribution is well described by three main 
components including a two component main cluster A2163-A with mass ratio 1:8, and its cluster 
satellite A2163-B. The bimodal mass distribution in A2163-A is similar to the galaxy density distri- 
bution, but appears as spatially segregated from the brightest X-ray emitting gas region. We discuss 
the possible origins of this gas-dark matter offset, and suggest the gas core of the A2163-A subcluster 
has been stripped away by ram pressure from its dark matter component. The survival of this gas 
core to the tidal forces exerted by the main cluster lets us infer a subcluster accretion with a non-zero 
impact parameter. Dominated by the most massive component of A2163-A, the mass distribution of 
A2163 is well described by a universal Navarro-Frenk- White profile as shown by a one-dimensional 
tangential shear analysis, while the singular-isothermal sphere profile is strongly ruled out. Comparing 
this cluster mass profile with profiles derived assuming intracluster medium hydrostatic equilibrium 
(H.E.) in two opposite regions of the cluster atmosphere has allowed us to confirm the prediction of a 
departure from H.E. in the eastern cluster side, presumably due to shock heating. Yielding a cluster 
mass estimate of M500 = 11.18^i['4g x 10^^h~^MQ, our mass profile confirms the exceptionally high 
mass of A2163, consistent with previous analyses relying on the cluster dynamical analysis and Yx 
mass proxy. 

Subject headings: cosmology: observations - dark matter - galaxies: clusters: general - gravitational 
lensing: weak - X-rays: galaxies: clusters: individual (A2163) 



1. INTRODUCTION 

Galaxy clusters are the largest self-gravitating systems 
in the universe, residing at the intersection of large- 
scale filamentary structures of the cosmic web. Ac- 
cording to the hierarchical structure formation scenario 
based on a cold dark matter (CDM) paradigm, mass ac- 
cretion flows onto clusters are still ongoing and high- 
mass ratio cluster mergers, the so-called major mergers, 
sometimes occur. Major cluster mergers are among the 
most energetic events in the universe, releasing amounts 
of gravitational energy as large as 10^^~^^erg. Spa- 
tially resolved X-ray observation has revealed to us that 
this collision energy partly dissipates in the intracluster 
medium (ICM) through shock heating and turbulence, 
yielding a complex ICM brightness and th ermal structure 
(see, e.g.. IMarkevitch fc VikhlininI 120071 and reference 
therein) . The dynamics of cluster collisions arc however 
dominated by the cluster mass distribution, which can- 
not be constrained from X-ray observations alone. Due 
to the coUisional nature of ICM, cluster mergers arc in- 
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deed expected to generate some transient decoupling be- 
tween spatial distributions of the CDM and hot X-ray 
emitting gas. Moreover, mass estimates relying on gas 
properties may be biased by merger-induced perturba- 
tions of the ICM hydrostatic equilibrium (H.E.) within 
individual colliding clusters. 

Weak-lensing distortions of background galaxy images 
provide us with a unique opportunity to reconstruct the 
distribution of matter in clusters without a ny assumption 
of ma s s model and dynamical states (e . g.. iKaiser et al. I 
T99I IBartelmann fc Schneideri [2OOII : ISchneideii 120061: 



Okabe fc Um etsul 120081: lOkabeet al.ll2010bD and mea- 
sure cluster masses (e.g., Grav et al.ll2002l: Gavazzi et al 



200I iBardeau et ai.ll2007t iHoekstral l2007i lOkabe et al 

2010bt lUmetsu et all 1200911 Therefor e, a weak 

lensin g anal ysis of merging cl usters (e.g.. iClowe et al] 

2OOI [Mahd avi et al. 200 71: lOkabe fc Umetsul 120081: 



Okabe et al.i i2010bi: iMerten et al.l I2011D is an obser 
vational breakthrough in measuring mass distribution, 
thereby providin g complementary ii iform ation to X-ray 
measurements. lOkabe fc Umetsul ()2008D conducted a 
systematic study of seven merging clusters, represent- 
ing various merging stages and conditions, based on 
a joint weak-lensing, optical photometric, and X-ray 
analysis, and revealed that the mass and optical light 
of member galaxies arc similarly distributed in merg- 
ing clusters regardless of their merging stages, but 
the mass distribution in merging clusters is highly ir- 
regular, and is quite different from the ICM distri- 
butions. It indicates that the gas and mass evolu- 
tions under cluster mergers are different. This fea- 
ture is confirmed by weak-lensing studies of 30 clus- 
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ters as a collaboration of "L ocal Cluster Substructure 
Survey" (lOkabe et al.l l2010bl). Thu s , a joint analysis 
of clu s ters (e.g., Mahdavi et alj 120081: lOkabe fc Umetsul 
200a lOkabe et all l2010d: iKawaharada et alJ 120101 : 



Umetsu et alJl2Q10l: IZtiang et alJl2010f ) will yield a com 



prehensive and quantitative understanding of cluster 
merger physics involved in the structure formation. 

Located at a redshift of 0.203, A2163 is a rich. X-ray 
luminous and hot galaxy cluster showing various signa- 
tures of ongoing merger events, including irregular opti- 
cal and X-ray morphologies (see, e.g. Elbaz et al. 1995 ; 



Markevitch fc Vikhlininll2001l:IMaurogordato et al.ll200S : 



Bourdin et al.l 120111) . and a pr ominent radio halo emis- 
sion (jFeretti et al.l l200ll I2004D . This cluster has been 
known to exhibit a clear spatial segregation between a 
bimodal galaxy dis tribution and a more centra lly peaked 
ICM morphology (jMaurogordato et al.l 120081 hereafter 
M08), while the projected distance separating its main 
mass and gas centroids has been measured as one of the 
largest in a sample of 38 clusters analyzed from both 
stron g galaxy lensing and X-ray imaging (jShan et al.l 
120101) . As revealed from spectroscopic and photometric 
analyses in M08, the galaxy distribution in A2163 can 
be separated into two components: the massive cluster 
A2163-A and its northern companion A2163-B, and the 
main cluster component itself showing a bimodal mor- 
phology with two brightest galaxies (BCGl and BCG2). 
In a more recent analysis o f Chandra and XMM-Newton 
data. JBourdin et al.l (J20111 hereafter Bll) evidenced the 
westward motion of a cool core across the E-W elon- 
gated atmosphere of the main cluster, A2163-A. Located 
close to the second galaxy over-density, this gas 'bul- 
let' appears to have been spatially separated from its 
galaxy component as a result of high-velocity accretion. 
From gas brightness and temperature profile analysis 
performed in two opposite regions of the main cluster, 
Bll further showed that the ICM has been adiabatically 
compressed behind this crossing "bullet" . 

Characterized by an exceptionally high ICM tem- 
perature first meas ured in X-ray (fc^T = 13.9 keV; 
lArnaud et al.l I1992[) then confirmed from its Sunyaev- 
Zeldovich (SZ) distortion (|Nord et al.ll2009[) . A2163 has 
been suggested to be exceptionall y massive from various 
X-ray and SZ analyses (see , e.g. IVikhlinin et al]|2009al : 
iPlanck Collaboration 1 120 lit) via the mass scaling rela- 
tion. Former weak-lensing analyses have been conducted 
using the V ery Large Teles c ope an d CF HT/Mega-Cam 
analys es of ICvpriano et al.l (I2004D and IRadovich et al.l 
(|2008[) . These analyses, however, did not exclude 
merger galaxies in the backgr ound shear catalog, as 
demonstrated by other studies (Broadhurst ct al."2005l: 
Okabc & Umetsu 2008; Umetsu & Broadhurst 200& 
Umetsu et al.l 120091 : lOkabe et al.l l2010bl: lUmetsu et aLl 
2O1O0. As shown bv lOkabe et al.l (|2010b( ). a contami- 
nation of member galaxies significantly dilutes lensing 
distortion signals, mainly in the central region, thereby 
yielding biased estimations on cluster parameters. 

In this paper, we conducted a weak-lensing analysis 
using Subaru/Suprime-Cam and made a secure selection 
to avoid a contamination of member galaxies in the shear 
catalog, combining CFHT/Mega-Cam data. The struc- 
ture of this paper is as follows: in Section[2]wc briefly de- 
scribe the weak-lensing analysis. The projected distribu- 
tions of mass, member galaxies and hot-gas are compared 



in Section [31 Section [J] is devoted to a one-dimensional 
shear analysis to measure cluster mass. In Sec [5l we 
perform, based on a two dimensional shear pattern, a 
multi-components fitting to measure three components 
revealed by M08. We compare weak lensing mass with 
dynamical and X-ray hydrostatic equilibrium masses in 
Section [51 In Section [71 we discuss quantitatively the 
physical process of making an offset between mass and 
gas distributions. Section |8] summarizes our results. The 
cosmological parameters of i^o = 100/i~^ Mpc"-"^ kms"'^, 
r^o = 0.3, and Aq — 0.7 are used in this paper. Given 
the cosmology, 1'. = 140.21 /i~^kpc. 

2. DATA ANALYSIS 

We retrieved Re image data from the Subaru archival 
data (SMOKA^). The data were reduced by the stan- 
dard imaging process of usi ng the reduction software 
for Suprim e-Cam, SDFRED (lOuchi et a l.| [200l. as de- 
scribed in lOkabe fc Umetsul (|2008f ) and lOkabe et al.l 
(j2010b[) . Astrometry calibration was conducted by fit- 
ting the final stacked image with the Two Micron All 
Sky Survey data catalog. The residual in astrometric fit- 
ting was lower than the CCD pixel size. The exposure 
time is 75 minutes and the seeing is 0'.'79; therefore, it is 
suitable for our weak-lensing study. Since no other band 
data are available, we use CFHT/Mega-cam g' imaging 
data to exclude unlcnsed galaxies in the shear source cat- 
alog via a color-magnitude plane. 

Our weak-lensing analysis is done using the IMCAT 
package provided by N. Kaiser (Kaiser et a l. 1995J ). 
We use the same pipeline as lOkabe et al.l (I2010airbl l 
with some modifications followed bv lErben et al.l ( 20011) . 
We first measure the image ellipticity, Bq,, from the 
weighted quadrupole moments of the surface brightness 
of each object and then correct the point-spread function 
(PSF) anisotropy by solving e^ = e„ - Pf^^(P;„,)^>''*, 
where Pa/s is the smear polarizability tensor and the 
asterisk denotes the stellar objects. We fit the stellar 
anisotropy kernel {Psm)al^^* ^^^^ ^^^ second-order bi- 
polynomials function in several subimages whose sizes 
are determined based on the typical coherent scale of 
the measured PSF anisotropy pattern. The median stel- 
lar ellipticities before and after the anisotropic correc- 
tion are (e^, e^) = (1.45 ± 0.02, 1.52 ± 0.01) x 10"^ and 
(ef"'*,e2''''*) = (-3.74 ±5.30, -5.07 ±4.73) x 10-^ re- 
spectively. Stellar ellipticities before and after the cor- 
rection and their pattern arc shown in Figures [1] and 
[21 respectively. We next estimate the reduced shear 
ga = 7a/(l - '^) = iPg)al^'i3 ^sing the pre-seeing shear 
polarizability tensor Pg, where we adopt the scalar value 
{Pg)a0 = Tr[Pg]5Q,^/2. We select background galaxies 
in the range of f^ ± (Jr* < rh < ?'h,max — 6.6 pixels, 
where r^ is the half-light radius, and f^ and Cr* are 
the median and standard error of stellar half-light radii, 
r/i* , corresponding to the half median width of circu- 
larized PSF. Then, we make a secure selection of back- 
ground galaxies in a color-magnitude plane in order to 
minimize a dilution of the weak-lensing signals caused 
by a contamination of unl ensed galaxies, main ly mem- 
ber galaxies. As shown in lOkabe et al.l (|2010b( ). the di- 

'' http://smoka.nao.ac.jp/indcx.jsp 

^ http://www.ifa. hawaii/kaiscr/IMCAT 



Subaru weak-lcnsing study of A2163: bimodal mass structure 



lution effect on the lensing signals is more pronounced 
at smaller radii because the ratio of the number den- 
sity of cluster galaxies to background ones rises towards 
the inner region. We use Subaru/Suprime-Cam Re and 
CFHT/Mega-cam g' images. The difference in filter sen- 
sitivity functions of the wavelength in the g' band be- 
tween these two instruments is negligible. Therefore, a 
combination with Subaru and CFHT enables us to effi- 
ciently exclude me mber galaxies inthe backgrou nd shear 
catalog. Following lUmetsu &: BroadhurstI ()2008[ ). we cal- 
culate the lensing signal as a function of color (Figure 
[3|), by averaging the tan gential distortion strengths (see 
also lOkabe et al.ll2010bD . We securely select background 
galaxiesin the range of 0.82 < {Rc~g') — {Rc — g')Rs < 2, 
-2 < (i?c - g') - {Re - g')RS < -0.38 and 22 ABmag < 
Re, where {Re — g')ns is the best-fit linear function as a 
magnitude Re of red-sequence galaxies; (Re — 5')rs = 
-0.056i?c + 2.812, as shown in Figure H Here, the 
magnitudes and colors are measured by MAG_AUTO 
(i.e., total magnitude) and MAG_APER (i.e., aperture 
magnitude) in SExtractor, respectively. By a selection 
in the color magnitude plane, the number density of 
background galaxies in the shear catalog decreases from 
~ 23 to ~ 12 arcmin"^. The mean redshift of these 
background ga laxies is estimated by matching the COS- 
MOS catalog (lllbert et all I2009D to calculate the aver- 
age lensing weight (Dk/Ds) = J^^dzdPwi,/dzDis/Ds, 
where Dg and Dis are the angular diameter distances 
between the observer and source (background galaxy) 
and lens and source, respectively, and dPwi,/dz is a 
probability function of redshift distribution. We obtain 
(As/A) ^ 0.7209. 

3. MAPS OF MASS, GALAXIES AND GAS 

According to the CDM paradigm, the coUisional ICM 
is expected to experience a spatial decoupling from galax- 
ies and dark matter during cluster collisions. Compar- 
ing hot X-ray emitting gas maps with the cluster mem- 
ber galaxy density and mass distributions has thus been 
known to testify to t he merging cluster dy namical states 
(see, e.g. iClowe eFal . 2006: Ok abe fc Umc tsu 200j). In 
order to investigate the collision scenario in A2163, we 
map out the cluster mass distribution and compare it 
with spatial distributions of the galaxy and hot gas com- 
ponents. 

As described in detail in lOkabe fc Umetsul ()2008D . 
we pixclizc the shear pattern into a regular grid using 
Gaussian smoothing kernel Wg{9) ex exp[— 0^/(?g] with 

9g = FWHM/V41n2. We adopt the smoothing scale of 
FWHM = 1.'33 due to the limitation of the number of 
background galaxies. We also use a statistical weight for 
each background galaxy in the context of 



1 



(1) 



s,« 



where ag^i is the rms error for shear measurement of zth 
galaxy and a is the softening constant variance. We set 
a = {Y.^l.il^f'^ - 0-44. The reduced shear at pixel 
position of 0„ is obtained by 



(ffa)(0. 



J2i'^g{(^n - di)UgS9a,i 
Y^i^gi^n - Oi)Ug^i 



(2) 



Then, we invert the pixelized reduced shear (Equation[5]) 
with a weight of the inverse of the variance at each pixel 
to the lensing convergenc e field, based on the Ka iser & 
Squires inversion method (jKaiser fc Squireslll993| ). 

We also map out distributions of optical luminosities 
and the number density of member galaxies with the 
same kernel with Ug = 1. The luminosity and density 
maps are sensitive to luminous and high-density struc- 
tures, respectively. We selected red-sequence galaxies 
with Re < 22ABmag as member galaxies. The absolute 
magnitudes for member galaxies are calculated from the 
apparent magnitude using the fc-correction for early-type 
galaxies. We assume that all cluster member galaxies in 
the catalog are at a single redshift. 

The resultant lensing k field is shown in Figures [5] and 
[6l with contours equi-spaced in units of la reconstruction 
error, Sk = 0.0404, above the Icr level. This mass distri- 
bution clearly exhibit a bimodal morphology in the cen- 
tral region (r ^ 3'.), with two peaks reaching significance 
levels of ^ lOtT and ~ 5.6cr, respectively. Interestingly, 
these two peaks coincide with the two galaxy overdensi- 
ties Al and A2 hosting the two brightest cluster galaxies 
revealed in M08, BCGl, and BCG2 (see also the top right 
and bottom left panels of Figure [6]). We hereafter refer 
to these primary and secondary mass peaks as MC and 
MW, respectively. The cluster mass distribution further 
exhibits some anisotropics at larger radii, mostly coincid- 
ing with peripheral optical clumps first revealed in M08 
and presently confirmed (see the bottom left panel of 
Figure [6|). The most significant of these coinciding mass 
and optical peaks is ~ 3.3cr at the optical clump, B. 

The mass distribution in the region of optical substruc- 
tures C, D, and E, discovered by M08, is similar to those 
of the number and luminous density, albeit low signifi- 
cance levels (Table [T]). We measure luminosities within 
2.' centering each k peak (Table [1]), where a background 
region of 20 — 22' centering BCGl is used. We find 
that luminosities and signal-to-noise ratios of k are corre- 
lated. We tried to measure model-independent mass for 
each optical subclump, following lOkabe et al. (2010a). 
However, since the result is sensitive to the choice of 
background region, we could not obtain reliable results. 
An over-density region is further revealed from the den- 
sity of galaxies whose colors are redder by ^ 0.6 at 
Re = 22 ABmag (see the bottom right panel of Figure 
[6|). The width of colors for these galaxies is ^ 0.15. Dis- 
tributed at the west of A2163, these optical structures are 
associated with mass clumps of ^ 2 — 3a. More generally 
speaking, the cluster mass distribution appears as spa- 
tially correlated with density distribution of its member 
galaxies, consistentl y with weak-lcnsing studies in other 
clusters (jOkabe fc U mctsu 2008; Okabe et al. 20101i). 

We next compare the cluster mass contours with X-ray 
surface brightness and temperature maps derived from 
XMM-Newton and Chandra data analysis in Bll. The 
left panel of Figure [7] shows us a wavelet denoised map 
of the overall cluster atmosphere, with mass contours su- 
perimposed. These two distributions clearly exhibit a 
spatial segregation, the X-ray emission appearing as uni- 
modal and E-W elongated, with an emission peak located 
between the two mass peaks. The right panel of Fig- 
ure [7] shows high-resolution details in a multi-resolution 
analysis of the Chandra image. This analysis reveals a 
secondary cluster X-ray core, XW, separated from the 
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main cluster emission peak, XC. As shown in Bll, this 
wedge-shaped feature (red-curved hne in the right panel 
of Figure [21) appears as a stripped cool core embedded in 
the hotter ICM of A2163-A and delimited in the west- 
ward direction by a cold front. The evidence of this cold 
front let us infer the westward motion of this gas 'bul- 
let' along the elongated atmosphere of A2163-A. Inter- 
estingly, we now observe that the crossing core appears 
as preceded by the secondary mass peak, from which it 
might have been spatially separated during a subcluster 
accretion onto the main cluster. This offset feature be- 
tween the X-ray cluinp and mass peak is coi isistent with 
other cold front clusters on a merging phase (iClowe et al.l 
[2006t lOkabe fc Unletsul [20081: lOkabe et al.l I2010bn . It 
suggests a possibility of ram pressure stripping of the 
accreted cluster core (see also Section [7|). 

It is further worth noticing that the main cluster mass 
peak, MC, also appears as slightly offset from the X- 
ray core, XC, so that the centroid of the gas core of 
XC and XW is at the intermediate position between the 
two mass peaks. This 'double offset' separating the gas 
and dark matter contents of both the infalling and main 
cluster cores presents some simi larities to the fas t ac- 
cretion observable in 1E0657-56 ([Clowe et al.ll2006[ ). the 
so-called "bullet-cluster". To our knowledge, no other 
bimodal merger has been known so far to exhibit such a 
significant offset between the gas and dark matter com- 
po nents of its two compoii ents (see, e.g. other examples 
in lOkabe fc Unlets^ [2008h . 

The northern subcluster A2163-B exhibits a spatial co- 
incidence between X-ray, mass, and member galaxies, 
consisten t with what is usually o bserved in pre-merging 
systems ([Okabe fc Umetsu[ [2008[ ) . Consistent with the 
lack of any interaction evidence found from ICM thermo- 
dynamics between A2163-A and A2163-B, this spatial co- 
incidence suggests that A2163-B is likely to be observed 
before interacting with A2163-A. 

4. TANGENTIAL DISTORTION ANALYSIS 

We conduct a tangential distortion study, which is a 
one-dimensional lensing analysis, in order to measure the 
total cluster mass. The tangential distortion component 
of the reduced shear, ga = (51,52), and the 45 degree 
rotated component for individual galaxies (ith galaxy) 
are obtained by 



9+,i = -gi,i cos 2ip - 52, i sin 2ip, 
5x ,i = -gi^t sin 2ip + g2s cos 2ip, 



(3) 



where Lp is the position angle in the counter clockwise 
direction from the first coordinate axis on the sky. Then, 
the profiles of 5+ and g^ are estimated with a statistical 
weight (Equation [TJ, as follows: 



(5a)(^„) = 






(4) 



Here, n denotes the nth radial bin 0„ with a given bin. 
The statistical error of 5^ in each radial bin is estimated 
as 



^) 



2 (E.«.)^ 



(5) 



where the prefactor 1/2 comes from the fact that ag^i in 
the rms is the sum of two distortion components. 



We fit the tangential distortion profile with the univer- 
sal profile proposed by Navarro, et al. (1996, hereafter 
NFW profile) and a singular isothermal sphere (SIS) halo 
model. The NFW halo mass is a prediction of numerical 
simulations based CDM model. The mass density pro- 
files over a wide range of masses are well described in the 
form of 

PNFw{r) = ' — -— ^, (6) 

where Ps is the central density parameter and Vg is the 
scale radius. The asymptotic inner and outer slopes for 
NFW mass density are p oc r^^ and r"'^, respectively. 
The three-dimensional mass profile within a radius, rg, 
at which the mean density is A times the critical mass 
density, pci-{z), at the cluster redshift, is expressed by 



with 



Mnfw(< ''a) = 3 to(ca). 



ra(x) =log(l + x) 



1 + a; 



(7) 
(8) 



The NFW halo mass is described by two parameters of 
the mass A/nfw(< '^a) and the halo concentration ca = 
?'a/?'s- 

The SIS halo model is a solution of the coUisionless 
Boltzmann equation, and is specified by one parameter, 
the one-dinicnsional velocity dispersion d^. as follows, 



Psis(?') 



1 



(9) 



27rG r2 ' 
A three-dimensional mass for the SIS model is given by 



Msis(<rA)--^rA. 



(10) 



The tangential distortion profile and the best-fit mod- 
els are shown in Figure [H We choose the central position 
determined by two-dimensional shear analysis which will 
be described in detail in Section [SJ We also measure the 
tangential shear profile with a center of BCGl and fit 
them, but the results do not change significantly. We 
can clearly find a curvature of the tangential shear pro- 
file. The curvature makes it difficult to fit the profile 
with the SIS model. The x^(d.o.f) for the SIS model 
is 27.85(7); therefore, the SIS model can be strongly 
rejected (5cr level) as a mass model. On the other 
hand, the NFW mass model well expresses the curva- 
ture of the profile. Indeed, x^ for the NFW model 
is 0.82(6). The best-fit NFW parameters are found in 
Table [2] The virial mass shows the massive cluster, 
Mvir = 24.251^;^^ x W^hr^MQ, where the virial over- 
density is Avir — 116.4. Our mass estimates are larger 
than those of a prev ious weak-lensing study (Table [5J 
[Radovich et al.l [2008D . using the shear catalog without 
excluding a contamination of member galaxies. 

5. TWO-DIMENSIONAL SHEAR ANALYSIS 

As shown in Figures [S] and [51 the projected mass 
distribution of A2163 is complex and two mass peaks 
are significantly detected in the central region. It is 
of prime importance for understanding cluster merger 
phenomena to measure masses of the main- and sub- 
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clusters. Furthermore , since numerical simula tions 
(jMeneghetti et all 120101: iBecker fc Kravtsovl [20lol ) and 



observations ( Okabe et al.ll2010bl) have shown that a tan- 
gential shear profile is affected by significant substruc- 
tures, taking into account substructures in modeling is 
also important for understanding such a lensing bias. In 
the mass measurement using the tangential shear pro- 
file (Section |4]), it is very difficult to distinguish which 
structure contributes in part to the tangential distortion 
signals, because the full lensing information from both 
the main and subclusters must be convolved to express 
the one dime nsional distortion p rofile with respect to a 
given center (lOkabe et al.ll201tel ). The two dimensional 
shear pattern, on the other hand, enables us to easily 
model lensing signals by a superposition of lensing sig- 
nals. In this section, we conduct two-dimensional shear 
fitting in order to measure the masses of three compo- 
nents (the sub and main components for A2163-A and 
A2163-B) revealed by M08 and Ell. 

We pixelize the shear pattern into a regular grid of 
If X If without any spatial smoothing procedure, whereas 
we adopted Gaussian smoothing in the map making (Sec- 
tion [3]) . The pixelized distortion signals and statistical 
weight, {ga){dn) and CTg(0„), in the nth pixel are esti- 
mated with a weight function Ui for each background 
source residing in the pixel (see also Equations, (jl]) and 
([Sp. The representative position for the nth pixel is also 
estimated with a weight function Ui. The x^ fitting is 
given by 



2 A^pixoi 
a, ,9=1 n 



{9p{0n)-9p{0n;pY'^°'"''^) 



(11) 



where p is the parameters and Cap is the error co- 
variance matrix of shape measurements in the form of 
CapiSn) = ^ap'^gi^n)- Here, S^p IS a Kronecker delta 
function and (T'^jOn) is t he statistical error of the pix- 
elized shear (jOguri et al.ll2010HWatanabe et al.ll2011f ). 

Wc first consider a single mass model of the NFW 
profile in order to compare the mass estimates by tan- 
gential shear measurement. We here treat the center of 
NFW mass {xc, Vc) as a parameter. In total, we use 
four parameters (M, c, Xc^ and j/c) for fitting. We adopt 
the Markov Chain Monte Carlo method with standard 
Metropolis-Hastings sampling. We restrict the sampling 
range of M^i, < 5 x IO^^/j-^Mq, Cvir < 20. We refer to 
the mean of the posterior probability distribution of each 
parameter. The resultant masses are consistent with the 
tangential shear measurements (Table [2]). The central 
position of mass is consistent with the peak position of 
the MC clump in the mass map. The MC clump, which is 
associated with BCGl, is therefore likely to be the main 
component. 

We next add a mass model for the mass clump MW to 
the main cluster. From X-ray and optical spectroscopic 
studies (M08 and Bll), the mass clump MW is likely 
to be a merging substructure in A2163-A. Since cluster 
substructure size is not determined by the virial theorem 
but by the strong tidal force of the main cluster (e.g., 
iTormeii et aI]|1998D we adopt the truncated SIS (TSIS) 
model (jOkabe et aLll2010af ) to describe the MW clump. 



The TSIS model is an extreme case of the truncation; 
mass density becomes zero at a radius r^. The TSIS 
mass profile is expressed as 



PTSis(^)=/Osis(f') 
= 



for r < rt 
for r > rt- 



(12) 



The subclump mass for the TSIS model is estimated as 



M 



(TSIS) 
sub 



2^ 
G 



rt- 



(13) 



The TSIS model is specified by one-dimensional velocity 
dispersion a^ and the truncation radius rt- We have an 
additional four parameters (cr„, rt and centers) for the 
TSIS model and a total of eight parameters for the fit- 
ting. We assume that the rcdshift of the MW clump 
is the same as that of the main cluster. Wc adopt 
\xc - a;poak| < 2.' and \yc - j/pcakl < 2.', where Xpoak and 
J/pcak are peak coordinates appearing in a weak-lcnsing 
mass map. The resultant mass and central positions are 
shown in Table [31 The centroid of the MW clump is 
consistent with the peak found in the mass map. The 
truncation radius, rt = 9.'7l2'9 '^ 1.4/i~^Mpc, is an in- 
termediate size compared to the virial radius of the main 
cluster Tvir = 16f6 ~ 2.3/i~^Mpc. Since the sub-cluster 
size after some impacts is significantly decreased by the 
strong tidal field, such a large truncation radius might 
suggest that the MW clump is a merging sub-cluster at 
the first impact. 

Next, we take into account the northern component 
A2163-B (M08 and Bll ), which contains luminous 
galaxies (optical clump B) and an X-ray emitting core. 
Wc consider two possibilities for the dynamic state of 
A2163-B; one scenario is the pre-merger phase that 
A2163-B is infalling toward A2163-A, and the other is 
that A2163-B has already undergone a merging event 
with the main cluster of A2163-A. Bll found no evidence 
of the close interaction between A2163-A and A2163-B, 
and therefore concluded that they are likely to be sep- 
arated more than rsoo aligning along the line of sight. 
However, in MOB, a filament of faint galaxies was de- 
tected along a north/south axis between A2163-A and 
A2163-B which might suggest a previous interaction be- 
tween the two components. Although this post-merger 
hypothesis is very unlikely from the X-ray approach, the 
origin of this faint galaxy filamentary structure is still an 
open issue. In this paper, we investigate two possibilities 
in the fitting. If A2163-B is physically separated from 
A2163-A, its mass profile is not likely to be affected by 
the strong tidal field of A2163-A and we therefore adopt 
the NFW model for the first scenario. Here, the halo 
concentration is ill constrained because the shear signals 
from A2163-B are smaller than those from A2163-A, and 
it is difficult in the environment of the massive cluster to 
find the curvature of the distortion profile, as shown in 
Figure [51 We therefore assume the mass - concentration 

relation of c^-„ = 7.85 (Mvir/2 x W^h-^M^y^'^^^ (1 + 
Zc)"°'^^ (|Duffv et al. l[20G8i) . Here, we assume a red- 
shift of A2163-B, Zc, is the same as that of A2163-A. 
This assumption is justified by both the dynamical anal- 
ysis (M08) and the ph otometric rcdshift estimates by 
iLa Barbera et all (|2004[ ). We also parameterize the cen- 
ter of A2163-B, using 13 parameters in total. The re- 
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sultant virial mass is A/vir3 = 2.43J;°;i5 ^ 10"/i~^Mq 

(Table O. The mass Msqo.b = ^■S'i.tofi x IO^/i-^'/q 
within a radius rsoo, at which the mean density is 500 
times the critical density, agrees well with X-ray esti- 
mates M5oo,s = (1.47 ±0.07) X lO^/i-^Afo (Bll). We 
changed the normalization of the mass -concentration 
relation by ±2 and conducted fittings, but the best fit 
of the virial mass for A2163-B changes only by -1-14% 
and —10%. They all are consistent within errors. We 
also found that the MC and MW masses are almost un- 
changed after adding the mass model of A2163-B. The 
resultant mass ratios of MC (main cluster of A2163-A), 
MW (sub cluster of A2163-A) and A2163-B arc found to 
be - 8 : 1 : 1. 

We consider the second scenario that A2163-B is an on- 
going merger with A2163-A. We here use the TSIS model 
for A2163-B, following the case of the substructure MW. 
The truncation radius r^ = 741.3^3428 kpc/i~^ is ob- 
tained, while the concentration of the NFW mass model 
is not well constrained. This is because the tangential 
distortion profile for the truncation model outside the 
truncation radius is proportional to the inverse square of 
the radius {g+ oc r~^), which is different from that of the 
NFW mass model for the main cluster. We obtain the 
substructure mas Mb = 1.44J:o:84 x lO^/i-^A/©. The 
mass discrepancy between the two models is small: 0.9a 
and 1.2cr with uncertainties of NFW and TSIS masses, 
respectively. The two mass models (NFW and TSIS) 
obtained solely by shear data are, therefore, acceptable 
for A2163-B. The total mass in each model is consistent 
within errors with tangential shear measurement, which 
indicates that a substructure effect is likely to be less 
significant. 

6. MASS COMPARISONS 

Due to its exceptionally high mass and dynamical ac- 
tivity, A2163 has long been an interesting test case for 
cluster mass measurements. We here compare the weak 



lensing cluster mass, M^ 



24.24 



-1-6.00 
-4.55 



X lO^^h-^Mp 



(see Section |4] and Figure H]), with dynamical and X-ray 
mass estimates provided by M08 and Bll. As previously 
mentioned, weak lensing mass measurements do not re- 
quire any assumption of cluster dynamical state. They 
consequently provide us a unique opportunity to investi- 
gate effects of the virial theorem hypothesis for member 
galaxies or the gas H.E. assumption on mass measure- 
ment precision. 

6.1. Mass estimate from galaxy dynam.ics 

We first compare weak-lensing mass with a dynami- 
cal one. MOB estimates dynamical total mass from the 
velocity dispersion, cri.o.s — 1400 kms~ , under assump- 
tions of the virial theorem, spherical symmetry, and no 
internal structure. The virial radius in dynamical mass 
estimation is defined as the three-dimensional radius of 
harmonic radius, corresponding to a coherent length of 
galaxy separations, which is different from that of the 
weak-lensing mass measurement. We here therefore com- 
pare M200 enclosed in the radius 7-200 within which the 
mean density is 200 times the critical density pcr- Dy- 
namical mass, M200 = (27.3 ± 2.8) x 10^'^MfD/t~\ ex - 
trapolated with a Hernquist profile (iHernguist I Il990l ). 
is higher than the weak-lensing mass measured by both 



one- and two- dimensional analyses (Table [5]) . The dis- 
crepancy might be due to a difference of models. We 
would require a more detailed dynamical model consider- 
ing u ncertainty of the anisot ropy in the velocity distribu- 
tion (jLokas fc Mamonll2001| ). On the other hand, we cal- 
culate one-dimensional velocity dispersion correspond- 
ing to the virial mass for the NFW model and obtain 

CTiD = T-^ii- i?(zc)A^(j. /2 ~ 1459 kms~^ assuming an 
isotropic velocity dispersion, which agrees well with the 
velocity dispersion along the line-of-sight (M08). They 
are higher than the one-dimensional velocity dispersion, 
n74.2Atfiit kms~\ of the SIS model which is ruled out 
at the 5(7 level by one-dimensional weak-lensing analysis. 
Wc note that this is a rough comparison, because the ob- 
served line-of-sight veloci ty dispersion is calculate d with 
a weight of mass density (|Binnev &: Mamonlll982|) . 

6.2. X-Ray m.ass estim.ates 
6.2.1. Mass estimate from ICM hydrostatic equilibrium 

We next compare X-ray masses using H.E. assump- 
tion. The impact of the ongoing cluster merger event in 
the central region of A2163 on hydrostatic equilibrium 
(H.E.) of the hot gas has been investigated in Bll. To 
do so, an average cluster mass profile has been extracted 
from density and temperature profiles of the overall clus- 
ter atmosphere, in addition to two supplementary profiles 
extracted assuming H.E. in the complementary regions 
located behind and ahead the gas core currently crossing 
A2163 (see also Sect. Eland Figure 7 of Bll). These pro- 
files have been compared with our NFW lensing profile in 
Figure [HI We first find that H.E. mass profiles extracted 
in the eastern cluster side and overall cluster strongly 
exceed the NFW profile at an overdensity radius of rsoo, 
and beyond. Interestingly, wc find instead that the H.E. 
mass profile extracted in the western cluster sector fully 
agrees with the weak- lensing profile in the radii range 
'"2500^'''500. Mass estimates corresponding to this analy- 
sis at the overdensity A = 500 are reported in Table [2j 
Adopting overdensity radii determined by each measure- 
ment method, we observe that the H.E. mass in all and 
eastern sectors are ~ 1.5 (^ 4cr level) and ~ 1.7 times 
(~ 5ct level) higher than the weak lensing mass, respec- 
tively, while the one in the western sector is in a good 
agreement with weak lensing one. From this compari- 
son, we may infer that H.E. assumption is only realistic 
in the western side of the cluster outer radii. The eastern 
cluster side is instead likely to have been shock heated 
by the ongoing subcluster accretion, yielding a strong 
departure of this cluster region and thus of the overall 
cluster atmosphere from H.E.. 

6.2.2. Mass estimate from the Yx proxy 

One concern in cluster cosmology is the search for X- 
ray mass proxies relying on well-calibrated scaling re- 
lations coup ling the cluster gas properties with total 
masses (e.g..lKravtsov et al.ll2006l:lVikhlinin et alll2009al: 



Arnaudetal.1 12 010': 'Zhang ct al.l 120081 : lOkabe et alJ 
2010c[ ). iKravtsov ct al. ( 2006 ) have proposed a new 
mass proxy, the so-called, quasi-integrated gas pressure, 
Yx = MgasT. This quantity has been suggested by 
numerical simulations to exhibit low intrinsic scatter 
regardless of the cluster dynamical state, in particu- 
lar since deviations of temperature and gas mass from 
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the normalization of the mass scahng relations partly 
anti-correlate during cluster mergers. The cluster mass, 
M500, has been estimated using the Yx proxy in Bll, 
assuming the M^qq-Yx scaling relation to be calibrated 
from hydrostatic mass estimates in a nearby cluster sam - 
plc observed with XMM- Newton (jArnaud et all l201Clf ). 
This mass estimate, -M500 = 13.09 



+ 1.40 
-1.68 



X lO^^Mn/i- 



yiclds a 20 % excess with respect to our weak-lensing 
mass. While marginally significant (l.lcr significance 
level), this mass discrepancy slightly exceeds the ^ 6% 
scatter of the Yx proxy s ugges te d by numeric a l simu- 
lations in iKravtsov et all ()2006[ ). lOkabe et all ()2010c| ) 
have investigated a dynamical dependence on mass ob- 
servable scaling relations, based o n XMM-Newton X- 
ray obse rvables (iZhang et al.l I2008D and weak- lensing 
masses (jOkabe et al.l l2010bf ). This study found the 
normalization of the M^qq-Yx scaling relation to be 
lower for disturbed clusters than for undisturbed clus- 
ters, contradicting the prediction of normalization in- 
dependence in cluster dynamical states. Mass esti- 
mates of A2163 relying on the undistur bed and dis- 
turbed cluster scaling relation calibrated bv lOkabe et al.l 

(12010^ yield M^qq'^''*^ ~ 14.53;^;^ x IO^Mq/i"! and 



M, 



(dist) 
500 



11.4 



+\i X lO^Mf. 



respectively. Inter- 



estingly, the mass estimate with the scaling relation for 
disturbed clusters fully agrees with weak-lensing mass, 
while the undisturbed cluster estimate yields a higher 
value, consistent with the Yx estimate of Bll. Since 
A2163 is an ongoing merger, the normalization for dis- 
turbed clusters gives a better result. I n order to con- 
strain the cosmological parameters (e.g, iVikhlinin et al.l 
l2009allH ). it is of prime importance to construct and cal- 
ibrate a low-scat ter mass proxy. R ecently, Millennium 
Gas Simulations (jStanek et al.l 120101 ) predicted that the 
temperature and gas-mass deviations are positive ly cor- 
related, which contradicts IKravtsov et al.1 ()2006l ) . It is 
also important to construct a low-scatter mass proxy 
using solely observational data, based on intrinsic co- 
variance measureme nt and principal component analysis 
(|Okabe et al.ll20103) . Further systematic studies, includ- 
ing both numerical simulation and observational study, 
are required for this. 

7. DISCUSSION 

From the inversion of a background galaxy shear pat- 
tern, we revealed a bimodal mass distribution in the cen- 
tral region of A2163 and various anisotropics at larger 
radii, the most significant of which corresponding to the 
northern subcluster A2163-B. Modeling the underlying 
three dimensional distribution of the cluster mass with 
three components further allowed us to constrain the 
mass ratio between the two central components to 1:8, 
and attribute comparable mass values to the central and 
northern subclusters. 

The spatial coincidence between X-ray peak, mass and 
member galaxy distribution in A2163-B suggests that 
A2163-B did not yet interact with A2163-A. While co- 
inciding with the bimodal galaxy distribution revealed 
in M08, the central mass distribution appears instead 
as spatially segregated from the brightest X-ray emit- 
ting gas region. In particular, the gas bullet suggested 
to cross the cluster atmosphere from the evidence of a 
cold front, in Bll, appears as preceded by the secondary 



mass peak to a projected distance oi ^ 2'. ^ 280/i~^ kpc 
along its westward direction of motion. Assuming an in- 
falling subcluster has experienced a transient separation 
of its gas and dark matter components as a result of ram 
pressure stripping, this configuration would require the 
gas core to have survived tidal distortions exerted by the 
main cluster, and the core temperature to be consistent 
with the subcluster mass. In the following, we discuss 
how these conditions allow us to put some additional 
constraints on the kinematics of the ongoing subcluster 
accretion. 

THE RAM PRESSURE STRIPPING CONDITION 

One possibility for explaining the offset between mass 
and gas distributions is the ram-pressure stripping: if the 
ram-pressure force on the X-ray core is stronger than the 
gravity around the core region, the gas is stripped from 
its gravitational potential. This condition is expressed 
as 



GM{< rcoic)Pc 



< A{Trr'^^^^psurv'^) I -7rr: 



3 

core 



(14) 



(jTakizawal I2006D . Here, r<.ore is the radius of the X-ray 
core, M{< rcore) is the spherical mass of the subcluster 
within the radius rcoro, Pcorc,sur are the density of the 
core and its surrounding gas, respectively, and v is the 
velocity of the core in the center-of-mass frame. A is 
a fudge factor at an order of unity. This is why the 
ram-pressure stripping is not effective due to a Kelvin- 
Hclmholtz instability, magnetic fields, and shock heating. 
Bll shows a lower limit of the density jump, Pcoro/Psur = 
1.28, at the west edge of the core (red-curved line in 
Figure III). 

We consider the ram-prcssurc stripping condition on 
the sub cluster. The substructure mass inside rcore is 
calculated from the TSIS model obtained by the weak- 
lensing analysis. Since Mtsis oc r, the condition (fH)) 
is free from the core size. Although we consider the in- 
ternal structure of gas, the resultant condition does not 
significantly change. Figure [10] plots the ram-pressure 
stripping condition Pram/-Pgrav — 1 as a function of in- 
fall velocity 1;, where Pram/-Pgrav gives the ratio of the 
right-hand side to the left-hand side of Equation ([T^ . 
We adopt A = 0.5 and 0.9. M08 found the gradient 
of line-of-sight velocity wios ^ 1250 km s~^, giving the 
lower-limit of infall velocity. The infall velocity is ex- 
pected from the masses to have an order of 



2G(Af„iai„ + M,ub) 



''"sub 



1/2 



2100 km s" 



(15) 



(jRicker fc Saraziiil 120011) . where M and r are the virial 
mass and radius, respectively. We here used observed 
truncated mass and radius for the sub cluster. The case 
of ^ = 0.9 satisfies the ram-pressure stripping condition 
for all velocities above Sv\os. The case oi A = 0.5 also 
satisfies the condition around 2100 km s~^. Therefore, 
the X-ray core initially associated with the sub cluster is 
easily stripped away from its central region. However, we 
must keep in mind that the estimated density jump might 
give the lower limit, due to the deprojection method we 
applied. In this case, a requirement of the ram pressure 
condition gives Pcore/Psm < 5(74/l)(ii/2100 km s~^)^. 
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CONSTRAINTS ON IMPACT PARAMETER FROM 
SURVIVAL CONDITION 

Wc next discuss the destruction process of the gas core 
by the stro ng tidal field of the main cluster. Numerical 
simulation (jTormen et al .112004 ) has shown that the life- 
time for a gas satellite is shorter than that for a dark 
matter one. In particular, the lifetime for a gas core de- 
coupled from dark matter potential is much shorter than 
that of a dark matter satellite. 

We roughly estimate the gas core disruption by the 
tidal force of the main cluster. As shown above, the 
gas core could easily be stripped away from the subclus- 
ter's core region. We therefore consider the stripped core 
which is composed of the gas, to be free from the gravity 
of the subcluster. The gas model for the core is adopted 
as a single power-law density profile as an approxima- 
tion. 



(16) 



ne{r) = 71e,0 — 



where the normalization ne,o and the slope a are the 
model parameters. We use the position of cold front 
ro = 0.'6 (red-curved line of the right panel of Figure [7]). 
The tidal radius at which the core density is truncated is 
estimated by the force equ ivalence between th e internal 
gravity and external tides (jTormen et al.|[T998() , 



n 



Ab 



Mcore(< n) 



-,1/3 



(2-ainM^ai„/ain6)Af„,ai„(6) 



(17) 



in the limit of r^ ^ 6 and Afcorc ^ A/main- Here rt is 
the tidal radius and b is the core position from the center 
of the main cluster; that is, an impact parameter. A is 
a fudge factor because hydrodynamic instabilities effec- 
tively dcstruct the gas core. We assume A = 1 in order 
to consider the tidal disruption only. We also assume 
that the density profile within a tidal radius dose not 
modify before and after the tidal stripping. The main 
cluster mass profile uses the best-fit NFW model in the 
case of NFW-f TSIS4-NFW (Tabled. We first calculate 
the minimum tidal radius rt^mm as a function of the den- 
sity normalization and its slope and impact parameter. 
The left panel of Figure [TT] shows a distribution of the 
minimum tidal radius of the gas core in the parameter 
plane of the density normalization and the slope. The 
contours denote the minimum tidal radius in the range 
of 0.'02 — O.'l. The point denotes the normalization and 
slope parameters for the core XW at the cold front. In 
the parameter plane, the minimum tidal radius is much 
smaller than the observed radius r ^ 0f6 ^ 120ft.yQ^ kpc, 
indicating that we cannot observe a remnant of the gas 
core if A2163-A is a system of close encounter. 

We next try to constrain the impact parameter, b, so 
as to realize the observed radial size of the core XW 
{rt = tq) through Equation. ([TT]). The right panel shows 
a distribution of the impact parameter in the parameter 
plane of the density normalization and the slope. The 
contours present the impact parameters (0.2 — 0.4)rvir to 
explain the core size. The observed X-ray core requires 
the impact parameter b ~ 0.25rvir ~ 560 kpc/i~^. If a 
subcluster collides into the main cluster with this impact 
parameter, the central region of the main cluster is sig- 
nificantly affected by cluster merger. It does not conflict 



with the disturbed core structures and the complex tem- 
perature distribution (Figure [T]). Since this estimation 
is based on the gravitational process only, we need to 
keep in mind a possibility that hydrodynamic instabili- 
ties shorten the lifetime of gas. We emphasize that this 
method demonstrates only one of the methods for joint 
X-ray and weak-lensing analysis. 

TEMPERATURE COMPARISON 

Assuming the cool gas core XW and the secondary 
mass peak MW represent the separated components of 
a formerly accreted subcluster, the intrinsic tempera- 
ture of XW should agree with the expectation of a 
self-similar subcluster temperature corresponding to the 
mass of MW. Given our multi-component mass distribu- 
tion (NFW+TSIS+NFW; see Table 2), we discuss the 
consistency of this hypothesis with projected ICM tem- 
peratures derived near the cool gas core XW in Bll. 

The cool gas core XW being located a relatively close 
projected distance from the main cluster emission peak 
XC, average ICM temperatures measured along a line of 
sight intercepting XW would consist of a linear combi- 
nation of the main and subcluster temperatures. These 
temperatures can be predicted within the radial range 
[0.2-0.5] rsoo, fro m the M500 - T[o.2-o 5] rsoo scaling rela- 
tion calibrated bv lOkabe et al.l (j2010cl ). in a sample of 12 
local clusters analyzed in both X-ray and weak-lensing. 
To perform this estimation, we first measure the main 
cluster mass, Msoo^main = 7.59ti:5? ^ lO^^h'^MQ, by 
muhi-component analysis (NFW-l-TSIS-f-NFW), within 
the overdensity radius rsoo. We further determine the 
subcluster mass, Msoo.sub, from its current mass, Mgub = 
Afmain/S, assuming the subcluster mass profile to have 
followed a universal NFW radial distribution before the 
ongoing accretion, and the profile concentration to be 
related t o M,^„,h by the halo mass-concentration depen- 
dence of iDuffv et al. I (|2008f ). These assumptions would 
yield a cluster and subcluster temperature prior to the 
accretion of fc^T = 10.lli:'lll;jtij keV and ksT = 

^■'''-0 9-0 3^1 2 kcV, respectively, where the first, second 
and third errors are from mass measurement error, nor- 
malization error of scaling relation, and intrinsic scatter 
in the relation. Assuming the main and subcluster emis- 
sivities to be comparable from a rough analysis of the 
X-ray image and weighting these temperatures following 
a scheme proposed in Mazzotta et al. (200 J), an average 
ICM temperature measured in the direction of the cool 
core would reach a value of fc^T ~ 5.7 keV. Despite be- 
ing consistent with uncertainties in our mass estimates, 
this value is considerably lower than projected temper- 
atures measured in the cool gas core XW in the XMM- 
Newton temperature map and the Chandra temperature 
profile of Bll (fc^T ~ 9 ± 1 keV). Among the possible 
origin for this discrepancy, the cool gas core might have 
been heated from its virial temperature while crossing 
the main cluster atmosphere, possibly due to mixing with 
the main cluster atmosphere or a reverse shock, presum- 
ably at an earlier stage of its accretion and prior to the 
formation of the cold front. 



In summary, the density of the main cluster atmo- 
sphere and the velocity of its presumably free-falling sub- 
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cluster are large enough to have separated the subcluster 
gas and dark matter components through ram-pressure 
stripping. Moreover, the survival of the subcluster gas 
core against tidal forces exerted by the main cluster im- 
plies the subcluster accretion to have occurred with a 
non-zero impact parameter. Comparing the average X- 
ray temperature in the direction of the cool core to its 
expectation from projection of the cluster and subcluster 
virial temperatures yields however a mild inconsistency, 
suggesting the cool core was partially heated while cross- 
ing the main cluster atmosphere. Deeper X-ray observa- 
tions should help better constrain the shape, tempera- 
ture and mass of the cool core and refine this accretion 



scenario. 



8. SUMMARY 



We presented the weak-lensing analysis of the merg- 
ing cluster A2163 using the Subaru/Suprime-Cam and 
CFHT/Mega-Cam data, measured cluster mass by one- 
dimensional tangential shear analysis, and measured 
three components' masses by multi-component analysis 
of the two-dimensional shear pattern. Based on comple- 
mentary X-ray, dynamical and weak-lcnsing datasets, we 
also discussed the centroid offset between weak-lensing 
mass and gas core. Our main results are summarized 
below. 

• The projected mass distribution shows a bimodal 
structure in the central part of A2163-A. The over- 
all mass distribution appears to be similar to the 
member galaxy one, whereas both mass and mem- 
ber galaxy distributions are completely different 
from the ICM one. This is consistent with a pre- 
vious study of seven mer ging clusters a t various 
dynamical state (jOkabc fc Umetsull2008f ). In par- 
ticular, we found a clear offset between the gas 
core associated with the cold front and sub-cluster 
mass peak. The offset was reported in all cold- 
front clu sters prev iously conducted by weak-lcnsing 
analvsis (iClowe et a l. 2006; Okabc & Umctsu 2003 
lOkabe et al.l l2010bf ). The gas core is also offset 
from th e main cluster ma ss peak, like the bullet 
cluster (jClowe et al.ll2006[ ). 

• A two dimensional shear analysis has enabled us 
to measure the mass of each of the three ma- 
jor components in A2163, including the two com- 
ponents of A2163-A, and the northern subclus- 
ter A2163-B. The central subcluster in A2163-A is 
well described by the TSIS model, giving the mass 
2.08tl^:^6 X lO"/i-iM0 ~ M„,ai„/8. 

• The mass of A2163-B is comparable to the sub 
cluster of A2163-A, which is in good agreement 
with the X-ray estimate of Bll, assuming H.E. The 
mass, member galaxies, and gas distributions are 
similar to one another, as r eport ed in pre-merging 
clusters fOkabe fc Umetsul ()2008[ ). suggesting that 
A2163-B did not yet interacted with A2163-A. 

• The gas bullet suggested to cross the cluster atmo- 
sphere from the evidence of a cold front, in Bll, 



appears to be preceded by the secondary mass peak 
to a projected distance of ^ 2.' ^ 280h~^ kpc. We 
show that the density of the main cluster atmo- 
sphere and the free fall velocity of an incoming sub- 
cluster with which mass is about one-eighth of main 
cluster virial mass are large enough to have sepa- 
rated the dark matter and gas component of this 
subcluster through ram-pressure stripping. Fol- 
lowing this scenario, the survival of the cool core 
against tidal forces exerted by the main cluster lets 
us infer that the subcluster must have been ac- 
creted with a non-zero impact parameter, reaching 
typical values of 6 ^ 0.25rvir ~ 560 kpch~^. As- 
suming the cool gas core and the secondary mass 
peak represent the separated components of the 
formerly accreted subcluster, the projected core 
temperature appears higher than expected from 
our estimates of the main and subcluster virial tem- 
peratures. Subsequent to its separation with dark 
matter, the crossing cool core may thus have been 
shock heated or disturbed by hydrodynamical in- 
stabilities yielding a partial mixing with the main 
cluster atmosphere. 

• Dominated by the most massive component of 
A2163-A, the overall mass distribution in A2163 
is well described by a universal NEW profile as 
shown by a tangential distortion analysis, while 
the SIS profile is strongly rejected (5(t confidence 
level). The virial mass for the NFW mass model is 



Mvir = 24.23l^;°g X lO^/^-ii 



• We compare the weak-lensing NFW mass with dy- 
namical and X-ray H.E. ones. The weak-lensing 
mass is lower than the dynamical one with assump- 
tions of the virial theorem and spherical symme- 
try mass distribution. This might be due to the 
dynamical mass model, because the linc-of-sight 
velocity dispersion expected from the NFW mass 
agrees well with the spectroscopic result (M08). 
The H.E. mass, M500, in the western sector show- 
ing no hot gas, is in good agreement with the weak- 
lensing one, whereas the one in the eastern sector 
showing hot gas is higher at the 5a level. This 
indicates that the merger shock heating leads to 
an overestimation of H.E. mass. The mass proxy, 
Yx , through the mass observable scaling relation 
for disturbed clusters (|Okabe et al.ll2010a ) gives a 
mass estimate similar to the weak-lensing mass. 
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Fig. 1. — Stellar ellipticity distributions before and after the PSF anisotropy correction. The left panel shows the raw ellipticity components 
(e*,e2) of stellar objects, and the right panel shows the residual ellipticity components (e^""' ,^2^' ) after the PSF anisotropy correction. 
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Fig. 2. — Distortion field of stellar cUipticities before (left) and after (right) the PSF anisotropy correction. The orientation of the sticks 
indicates the position angle of the major axis of stellar ellipticity, whereas the length is proportional to the modulus of stellar ellipticity. A 
stick with the length of 10% ellipticity is shown above the left panel. 
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Fig. 3. — Mean distortion strength (((?+, x)) over the radii of 1' < r < 18', as a function of color. The background samples are defined 
with galaxies redder or bluer than the red-sequence by at least the color offsets given by the x-label. The distortion strength is changed 
due to dilution by cluster members. The red circle and blue square points denote the mean tangential distortion {(g+)) for redder and 
bluer background samples, respectively. The green circle and square points denote the mean tangential distortion for the 45 degree rotated 
component ((^x )) for redder and bluer samples, respectively. The two vertical solid lines d enote our choices of the color cuts used to define 
the red /blue background galaxy samples shown in Figure U] The details are described in lUmetsu fc Broadhurs3 II2008I ') and lOkabe et ahl 
1201 OH') . 



12 



Okabe et al. 




IB 30 22 

R^ [AB mag] 



Fig. 4. — Color-magnitude diagram. Green points are member galaxies of which spatial distribution is shown in the top right and bottom 
left panels of Figurc[6| Two green dashed lines represent the width of the red sequence. The red and blue points are the background galaxy 
samples, used for the Icnsing distortion analysis, whose colors arc redder and blue than the red sequence, respectively. The magenta points 
denote a sample of galaxies in the right bottom panel of Figure \E\ 
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Fig. 5. — Subaru _Rc-band image (23.' X 23.') overlaid with contours of the K-ficld reconstructed from the weak-shear field, in units of Icr 
reconstruction error Sk = 0.0404. 
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Fig. 6. — Top left: reconstructed lensing K-field in units of la reconstruction error. The bimodal structure in the central region is clearly 
found. The significance levels of bimodal peaks are ~ lOcr and ~ 5.6cr, which are referred to as MC and MW, respectively. Top right : 
luminosity map for red-sequence galaxies in the Re band. Overlaid are the contours of the k field. Two crosses denote BCG positions. 
Bottom left : densi ty map for red-sequence ga laxies in units of arcmin"'^. The four density clumps other than the bimodal structure, 
discovered by M08 I IMaurogordato et al.l 120081) . are confirmed. All labels for the density clump are quoted from MOB. Bottom right : 
density map for galaxies whose color is redder by ~ 0.6 at Re = 22 ABmag than those of member galaxies. 
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Fig. 7. — Left: adaptivcly-snioothcd Chandra X-ray image (14.' X 14.') in the 0.5 — 2.5 keV band. Overlaid are the contours of the lensing 
K- field in units of la reconstruction error (Figure |6ll. Middle : temperature map from XMM-Newton data with the same mass contours 
as those in the left panel. Right: high resolution residual Chandra image obtained after subtracting large scales of the denoised wavelet 
transform of the X-ray image. Two components (XC and XW) are clearly found in the X-ray emitting core. Two crosses denote BCG 
positions. Red (color version) or dark gray (black and white version) curved line represents the cold front (CF; see Bll). Overlaid are the 
same mass contours above the Scr significance. The offsets between the X-ray core and two mass peaks are clearly found. 



TABLE 1 
Optical subclumps 


Names 
(1) 


iRj<2.') 
(2) 


S/N 
(3) 


C 
D 

E 


1.47 
0.94 
2.38 


2.0 

1.7 
2.7 



Note. — Column (1): names of optical subclumps identified by M08. Column (2): 
Column (3): signal-to- noise ratios appearing in the mass map (Figures [5] and |6j. 



the luminosity within 2f in units of W^^ LqH 



TABLE 2 

Best-fit Mass Profile Parameters 



Method 


(IQI^/i-IMq) 


Cvir 


(iQi^h-iMe) 


C200 


M500 
(lOi^ft-iJUe) 


C500 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


ID WL 
2D WL 


24.24+^-^^t' 
26.72tt:?9 


2 84+"-''-^ 
9 co+0.48 


20.32±igi 


9 1 «+('■''>''■ 

^■1^-0.50 

t-yy_o.49 


u.isl-ll 

ii.i3iliJ 


1 oc+u.yy 

1.00_Q 33 
1 07+0.20 

^^■•^'-0.30 


Dynamics (M08) 


- 


- 


27.3 ±2.8 


- 


- 


- 


X-ray (Bll) : H.E. 


- 


- 


- 


- 


n.29tlil 


- 


X-ray : H.E. west 


- 


- 


- 


- 


10 1^+0-42 

iu.io ; Q 21 
18.55_Q;g4 


- 


X-ray : H.E. east 


- 


- 


- 


- 


- 


WL (Radovich et al. 2008) 


15.4 ±2.8 


- 


12.6 ±2.1 


- 


8.4 ±1.4 


- 



Note. — The best-fit NFW parameters obtained by fitting to the one-dimen sional tangential distortio n profile (ID) and the two- 
dime nsional shear pattern (2D). For com parison, the masses de rived from dynamical IIMaurogordato et al.l2003 . M08), X-ray IIBourdin et al.l 
120111 . Bll), and weak-lensing analyses JRadovich et al.ll2008l ) are quoted. H.E., H.E. west and H.E. east are hydrostatic masses for all, 
western, and eastern sectors, respectively. X-ray mass estimations uses rsoo determined solely by X-ray analysis. Weak-lensing masses 
are obtained using the shear catalog without excluding a contamination of member galaxies. Column (1): method. Column (2) and (3): 
virial mass in units of IO^^/i^^Mq and the concentration parameter for the NFW model. Column (4) and (5): M200 for the overdensity 
of A = 200 and concentration parameter. Column (6) and (7): M500 for the overdensity of A = 500 and concentration parameter. 
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Fig. 8. — Profiles of the tangential shear component (top panel), g_|_, and the 45 degree rotated component (bottom panel), gx- The 
solid and dashed lines are the best-fit NFW and SIS models, respectively. The SIS model is strongly disfavored (~ 5a level) for a cluster 
mass profile. 



TABLE 3 

Multi-component Analysis of the Two-dimensional Shear Pattern 



Parameters 


NFW-I-TSIS 


NFW-fTSIS+NFW 


NFW-fTSIS+TSIS 


(1) 


(2) 


(3) 


(4) 


-^vir,inaiii 


20.9Sttil X W^^h-^MQ 


16.46+^-^;; X IO^/i-^Mq 


18.79+e-^g X lo"/i-^M0 


'^vir,main 


r, oa+0-70 


o 44+0.78 


, 14+0.74 


Mmw 


2.13+03^ X lO^^h-^MQ 


2.08+S 8? X lO"?i-iM0 


2.09+° ;J^ X lo"/i-iM0 


'"t,MW 


1.39l°1>-iMpc 


1.35+St>-'Mpc 


1.37+Sf3/i-iMpc 


Mb 


- 


2.43+°';? X lO^^h-HlQ 


lA4t°oit X W'^h-'M^ 


(": (5)0, main 


{243.95stVoil -e.urtlZl) 


(243.958;°:r3. -6.14612:^3) 


(243.958«:™t -6.147l°r3) 


(a, 5)c,MW 


(243.905inS5. -6.125tS:r5) 


(243.9051°;™. -6.125tS:ra) 


(243.9051°:°™, -6.125l°;ra) 


[a, 5)c,B 


- 


(243.9691°:™, -6.042lS:ro5) 


(243.968l°:°„S^, -6.042lS:SS|) 



Note. — The best-fit parameters obtained by fitting the two-dimensional shear pattern with a multi-component model of clusters: 
NFW parameters (Mvir main, ^vir main) for the main cluster (MC) at the virial over-density, the truncated-SIS (TSIS) mass (Mmw) for 
the subcluster (MW), its truncation radius (rj mw) iii units of Mpch~^ , and the mass of A2163-B (Mb). The best-fit halo centers are in 
units of degree. Column (1): parameters. Column (2): NFW and TSIS models for the main cluster (MC) and for the substructure (MW), 
respectively. Column (3): NFW, TSIS and NFW models for the main cluster (MC), the subcluster (MW) and A2163-B, respectively. 
The mass Mb is the virial mass. Column (4): the NFW model for the main cluster (MC), the TSIS model for the subcluster (MW) and 
A2163-B, respectively. Mb is the mass within the truncation radius. 
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Fig. 9. — Comparison of NFW lensing mass and hydrostatic masses. Black color rep resents the NFW len sing mass profile. Green, red 
and blue colors represent the hydrostatic mass in overall, east and west sectors (Bll; IBourdin et al.ir2011l) . respectively. The solid and 
dashed lines are the best fit values and the 68% CL uncertainty errors, respectively. The arrows denote the ovordcnsity radii r25oo and 
rsoo determined by weak-lcnsing analysis. We use ^70 = Hq/TO km s~^ Mpc"-'^ = 1.0. 
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Fig. 10. — Ram pressure stripping condition of the subclustor core. A is a fudge factor associated with hydrodynamical instabilities or 
gas heating^ The vertical line denotes the gradient of the line-of-sight velocity (M08). The arrow denotes an order of the infall velocity 
(Equation IllSt ') expected from masses for main and sub clusters, measured by two-dimensional weak-lensing analysis. 
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Fig. 11. — Left ; ininimum tidal radius in the parameter plane of core density normalization, ne(ro) = riefi at the cold front r = ro 
and slope a. The contours represent tidally truncated core sizes of 0.'02,0f05, and Ofl in units of arcmin. The point denotes the observed 
parameters for the core XW. Right: impact parameter to realize the observed core size. The contours are (0.2 — 0.4)rvir. The point is the 
same as that in the left panel. 



